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Impulsive pp-waves are commonly described either by a distributional spacetime metric or, alterna- 
" tively, by a continuous one. The transformation T relating these forms clearly has to be discontinuous, 

■ which causes two basic problems: First, it changes the manifold structure and second, the pullback of 

' the distributional form of the metric under T is not well defined within classical distribution theory. 

Nevertheless, from a physical point of view both pictures are equivalent. In this work, after calcu- 
» I 1 lating T as well as the "Rosen" -form of the metric in the general case of a pp-wave with arbitrary 

O l' wave profile we give a precise meaning to the term "physicially equivalent" by interpreting T as 

, the distributional limit of a suitably regularized sequence of diffeomorphisms. Moreover, it is shown 

that T provides an example of a generalized coordinate transformation in the sense of Colombeau's 
generalized functions. 
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I. INTRODUCTION 



o 

7— I 

00 
On 

^ , Plane fronted gravitational waves with parallel rays (pp-waves) are spacetimes characterized by the existence of a 
I covariantly constant null vector field, which allows to write the metric in the form [1] 
5— ( . 

ds 2 = Hdu 2 + 2d(d( - 2dudv , (1) 

J> '. - 
. \ where H(u, £, £) is a function of the the retarded time coordinate u and the transverse complex coordinates Q and Q. 

■ Impulsive pp-waves may now be defined by setting [2] 

H(u,CX) = /(CO <*(«)> (2) 

where / denotes a smooth function of £, £ and S is the Dirac-(5-distribution. Such spacetimes arise in a very natural 
manner as the ultrarelativistic limit of Kerr-Newman black holes [3-5] and multipole solutions of the Weyl family as 
recently shown by Podolsky and Griffiths [6] . Moreover these geometries play an important role in the description of 
scattering processes at the Planck scale [7,8]. 

While the form (1) of the metric very clearly demonstrates the nature of the impulsive wave, i.e. that the space- 
time is flat everywhere except for the null hyperplane u — where the 5-like shock is located, it has the obvious 
disadvantage of involving distributional coefficients. However, due to the simple form of the metric the curvature 
tensor can be calculated and is itself proportional to the ^-distribution. Moreover, the geometry of impulsive pp-waves 
was described recently [9,10] by the authors entirely in the distributional picture. The key idea was to use a careful 
regularization procedure which allowed to handle the nonlinear singular geodesic and geodesic deviation equations in 
a mathematically satisfactory way. 

On the other hand, impulsive pp-waves are frequently described by a different spacetime metric which is actually 
continuous (see [2,11,12] and, for the general case, [13]). In the special case of an impulsive plane (linearly polarized) 
wave, i.e. / = 1/2 (C 2 + C 2 )i h takes the form 
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ds 2 = (1 + u+fdX 2 + (1 - u+fdY 2 - 2dudV . (3) 

where u + denotes the "kink" function u + = 9(u) u. This form of the metric has the advantage that only the curvature 
tensor involves distributions while the metric can be treated "classically." Moreover, if one constructs impulsive 
pp-waves by glueing together two copies of Minkowski space along the hypersurface u = with a coordinate shift 
(according to Penrose's scissors and paste approach [2]) a theorem of Clarke and Dray [14] ensures the existence of a 
C^-atlas in which the metric components are continuous. 

Clearly a transformation relating the metrics (1) and (3) cannot even be continuous, hence strictly speaking it changes 
(even the topological structure of) the manifold. In the case of a plane wave this discontinuous change of variables is 
given by the so called Penrose junction conditions [2] (for the general vacuum case again see [11]) 

V2( = X(l + u+) + iY(l - u+) 

v = V + \x 2 (u + + 6{u)) + l -Y 2 {u + 9{u)) (4) 
u = u. 

However, the two (mathematically) distinct spacetimes are equivalent from a physical point of view, i.e. the geodesies 
and the particle motion agree [15]. More precisely, if one computes the geodesies for the continuous metric and then 
(formally) applies the discontinuous transformation one obtains exactly the geodesies as derived in the distributional 
picture in [9]. 

Nevertheless, it has to be noted that the actual calculation of the continuous metric (3) from the distributional form 
(1) as well as the transformation of the geodesies involve undefined nonlinear operations with distributions which are 
dealt with by using ad hoc multiplication rules; in this case 9 2 = 9, 9u+ = u+, u+5 = = u 2 + 5 and 95 = 1/2 6. Such 
a seemingly harmless operational approach to introducing nonlinearities into linear distribution theory may work out 
in certain cases but in just as many instances it will lead seriously astray. Using the above, i.e. exactly the "rules" 
appearing in the situation under consideration we get {99)8 = 95 = 1/2 5 =/= 1/45 = 9(95). For examples of the 
difficulties that may result from such manipualtions both in mathematics and in physics, cf. [16] and [17]. Especially 
in the context of nonlinear differential equations involving generalized functions (in this case: the geodesic equations) 
it turns out that the method of choice for determining which multiplications work out and which don't is to first 
regularize the singularities, then carry out the nonlinear operations and, finally, compute the distributional limits, if 
they exist (cf. [18], p. 180). 

By following this route, in this paper we give a precise meaning to the above-mentioned "physical equivalence," inter- 
preting the discontinuous change of variables as the (distributional) limit of a family of smoothened transformations. 
The key fact is to note that the coordinate lines in the new variables defined in (4) are exactly the geodesies of the 
metric (1) with vanishing initial speed in the (— , £— and v— directions. 

More precisely, we start with the distributional form (1) of the metric, which we turn into a sandwich wave by 
regularizing the 5 distribution in a completely general manner. Then we carry out a transformation analogous to 
(4) and calculate the distributional limit of the metric to arrive at the continuous form, which in the special case 
considered above precisely agrees with (3). In calculating this limit, which is totally independent of the form of the 
rcgularization (hence "natural"), we make use of the results in [9] and [10]. In the sandwich wave - picture both forms 
of the impulsive wave, i.e. the metrics (1) and (3), arise as (distributional) limits in different coordinate systems (see 
also the remarks in Sec. V). 

II. THE SMOOTHENED TRANSFORMATION 

Since our considerations depend heavily upon the results of [9,10] we shall stick to the real coordinates used there 
and start with an impulsive pp-metric of the form 

2 

ds 2 = f(x l ) 5(u) du 2 - dudv + ^(cfa i ) 2 , (5) 

i=i 

where x l = (x, y) denotes the transverse coordinates spanning the wave surface. The first step is now to regularize the 
^-distribution by a sequence of (smooth) functions p £ (0 < e < 1). Since we are interested in the most general result 
we will only assume the following minimal set of conditions on the regularization (so called generalized 5-function) 



(a) SUpp(p £ ) C [-£,£] , 
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(b) / p e (x)dx^l (s — > 0) and 



dx < C for small e . 



(6) 



(C) J \Pe{x) 



Note that condition (a) is choosen in order to avoid technicalities in the calculation of the distributional limits which, 
however, remain valid if (a) is replaced by 



(a') supp(p e ) - {0} (e - 0) . 
Hence we are left with the following sequence of sandwich waves 



2 



J\2 



ds £ = f(x l ) p £ (u) du - dudv +> (dx 1 ) 



»=i 



(7) 



(8) 



The corresponding geodesic equations are given by 

1 



xl(u) = -d i f(xi(u))p £ (u) 
v £ (u) = f(x{{u))p £ {u) + 2dif(xi(u)) x l £ (u)p £ (u), 



(9) 



where we have used u as an affine parameter. Since we are only interested in geodesies with the special initial 
conditions x l £ (— 1) = x l , v £ (— 1) = vq and x l £ (— 1) = = v £ (— 1), we adopt the notation xI(xq,u) and v £ (vo,Xq,u) f° r 
these geodesies with vanishing initial speed. They obey the following (implicit) set of equations 



x\ (Xq, u) — Xq + — 



dif{xi(xl,r))p £ (r)drds 



U US 

v e (v ,xi,u) =v + J f(xl(x^s))p E (s)ds + J J dif(xi{x%,r))x E {x%,r))p £ (r) dr ds 



(10) 



In [10] it was shown that these equations are uniquely solvable (even) within the framework of Colombeau's algebras of 
generalized functions [18-20], which provide a suitable setting for singular differential equations. Moreover the unique 
solution possesses the following distributional limit, independent of the concrete shape of the regularization 



xi{xl,u) -» 4 + -d i f(x 3 )u + 
v e (v , 3p Q , u)->v + f(xi) 9{u) + ^difixi) d\f(x 3 a ) u + . 

Let us now consider the transformation T £ : (u, v, x l ) — > (u, V, X 1 ) (where X % 
depending on the regularization parameter e according to 

x l = x\(X j ,u) 
v = v £ (V,X j ,u) 



(11) 

(X,Y)) given implicitly by (10) and 

(12) 



which is precisely analogous (and in the limit e — > and the special case of a plane wave actually reduces) to (4). 
We claim that T £ yields a coordinate transformation such that the new coordinates (u, V, X, Y) are constant along the 
geodesies given by (10). While the latter property follows directly from the construction, we still have to verify that 
T £ is a diffeomorphism in a spacetime region containing the shock-hypersurface u = 0. To establish this property we 
employ a global univalcnce theorem by Gale and Nikaido ( [21], Thm. 4) stating that any differentiable F : fl — > R n , 
where is a closed rectangular region in R™ is univalent (injective) if all principal minors of its Jacobian J(x) are 
positive. Since 



d( u, x 1 , x 2 , v ) 
d(u,X\X 2 , V) 



1 






dx x E 






du 


ax 1 " 




dxl 


dxl 


dxl 


du 


ax 1 " 


dX? 


dv e 




dv £ 


du 


ax 1 " 


ax 3 



(13) 
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we have to find estimates for 

dxl 



U S 

= 6i j+2j J{dmdd){x k e {X l ,r))^-{X l ,r)p £ {r)drd S . (14) 



If X 1 varies in a compact region K of R 2 and — l<p<u<ooit follows from a straightforward modification 
of the appendix of [9] that for small e the terms x £ (X l ,p) remain bounded, independently of X 1 , p and e. Let 

g(u) := sup{£ \§£j{X\p)\ : X k e K, -1 < p < «}. Then (14) gives 

i=l 



u 

\g(u)\<C 1 + C 2 J \g(s)\d s . 



so Gronwall's lemma implies the same boundedness property for ^jfj- Using these estimates it follows from (14) that 

for small e and for small u > (depending on \\didjf \\oo in a compact region and on C from (6c)) will remain 
arbitrarily small (for % ^ j) or arbitrarily close to 1 (for i = j), respectively. Hence all principal minors of (13) are 
indeed positive in a suitable rectangular region (independent of e) containing the shock hypcrsurface u — which 
establishes our claim. We note that this result in particular implies that for small e the geodesies (10) do not cross in 
this region. 

In the new coordinates (u, V, X, Y) the metric tensor takes the form 

2 2 

ds 2 e = -dudV + (2^(4^-4) - d j v £ )dudX j + Y,(djxidX j f , (15) 

i=l i=l 

where ' and di denote derivatives with respect to u and X 1 respectively and we have omitted the straightforward 
computation showing that g uu = 0. 

III. DISTRIBUTIONAL LIMIT 

Our next task is to compute the distributional limit of the metric (15) as e — > 0. To reduce the notational complexity 
we are going to suppress the explicit dependence of test functions on the variables X 1 and V, i.e. we will write <p(u) for 
ip e 2?(R 4 ) and will simply drop the integrations with respect to X % and V. This abuse of notation is admissible due 
to the uniform boundedness properties (in X % and V) established above and serves to reduce the number of integrals 
in the sequel by 3. We start out with the coefficient 

2 u us 

9ux =E(/ dtfixKX^s^pMds d.iX 1 + l -JJ dJ(xi(X\r)p £ (r)drds)) 

% = 1 -6 -6 -6 

u US 

-di(Y + J f(xi(X k ,r))p E dr + J J a i /(4(X fe ,r)4(^',r)p £ (r)rfrrfs) 

— e —e —e 

The first term of this expression can be written as 

2 u u us 

J2(Slj d l f(xl(X\s))p e (r)dr+ l - J dJ(xi(X\r))p e (r)dr d,J J cW £ (A fe , r)p £ {r) drds) 



(16) 



-^dif(Xi)6 —A 

For any test function ip we have 



oo u us 

(2A,p) = J p(u) J dJ(xi(X k ,r))p £ (r)drd 1 J J d l ]{x\{X k ,r)p £ (r) drdsdu 
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Splitting this integral into a sum of the form 

e u us oo u e s oo u us 

II- (17) 



-€ —S —S — € € — € — € —S € —6 5 —€ 



the boundedness arguments following (14) imply that the first two summands converge to 0. Also, by (11) in [9] and 
by (10) we obtain 

lim sup \di(f(xi{X k ,r))) - 3J(A fe )| = (18) 

|r|<e 

uniformly for X k in compact sets. For later use we note that by the same reasoning also 



lim sup \d i {d J f(xi{X k , r))) - d l0 f{X k )\ = (19) 
uniformly for X k in compact sets. Hence by a direct estimation the limit of the remaining term is 

2 

(d i f(Xi)d 1 dif(Xi)u+,ip), so the distributional limit of the first term in (16) is d 1 f(Xi)6+± J2 d i f{X^d i d 1 f(X^)u + . 

i=l 

Turning now to the second term in (16) we obtain from (18) and (6b): 

u 

d 1 J f{4{X\r))p £ {r) dr - 

— e 

Inserting (10) it remains to calculate the distributional limit of 

us r 

\J J di(dif(xi(X k ,r))p s (r) J d i f(xi(X k ,q)p s (q)dqdrds 

— e —e —e 
us r 

+ \j j dif(xi(X k ,r)p £ (r) I d 1 (d l f(xl(X k ,q))p E (q)dqdrd S . 
—e —e —e 

A splitting scheme as in (17) for both these terms, (18) and the fact that 

us r 

I j Pe ^ j ' P^)dqdrds -> 

e —e —e 

imply that each summand in (20) converges distributionally to \d l f{X^)d u f{X^)u+. 
An analogous argument holds for g u y ■ Summing up, we obtain 



(20) 



If we write 



then 



X)' d' 
g uX — ► and g uY — > for e — > . (21) 



2 u s 2 

gxx = J2(di(X i + \ j jd i f(xi(X k ,r))p E (r)drds)) 2 =:J2R i (22) 



US US 

Ri = Sj + J J d 1 (d 1 f(xi(X k ,r)))p e (r)drds+\(J J d 1 (d i f(xi(X k ,r)))p £ (r) dr daf 



—e —e 
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and by a similar argument as above it follows that the distributional limit of the derivative d u B is \{d\dif{Xi)) 2 u + . 
Since taking primitives (i.e. convoluting with 9) is separately continuous on the convolution algebra of distributions 
supported in an acute cone and since taking tensor products of distributions is separately continuous as well we get 

1 ?/ 2 
Ri ^ 5} + 9 n /(X>+ + -(cW^)) 2 ^ . 

Thus 

9xx ^ (1 + ^ n /(X> + ) 2 + l(a 12 /(^)) 2 4 (23) 
9yy (1 + \d 22 f{X=)u + ) 2 + \{d 12 f{Xi)fu\ (24) 

2 2 

Finally, we turn to g XY = 2Y J {d\x\{Xj))(d 2 x\{Xj)) =: 2 £ Inserting from (10) we have 

i=i i=i 



+ 



11 s US 

Sl = \j J d ^fH( xk > r )))Pe(r)drd.s+± J J d 1 (d 2 f(xi(X k ,r)))p e (r)drds 

— £ —e —e —e 

^ v ✓ \ v ■* 

^dl2 f(X0) U+ -.912 f(X3) U+ 

us us 

\J j d 2 {d i f{4{X k ,r))) Ps {r)drds j j d l (d i f(x> l {X k , r)))p £ {r) dr ds 

— e —e —e —e 

" v ' 

=:D 

We claim that 

D - d u f{Xi)d 2i f{Xi)u 2 + . (25) 

To establish this, let ip £ V and consider 



oo / u s 



(D,<p)-d u f(Xi)d 2i f(Xi) I u 2 y{u)du = j 1 1 1 d 1 (d i f(xl(X k ,r)))p £ (r)drds 

—e \-e -e J 

)oc 
<p(u) du- j d u f{X r )d 2l f{X3)u\{u) du 
o 

In order to show that this goes to 0, by a splitting similar to (17) and by the boundedness properties already established 
it suffices to prove 

oo /us 

d 1 (d i f(xi(X k ,r)))p e (r)drds]\ I I d 2 {d t f{x{{X k ,r)))p E (r) dr ds I <p(u) du - 



OO 

I ftj/^'jfe/^JttVu) du -> , 



which in turn is a consequence of (19) and (6). Therefore, 

5xy ^ (dnf(X j )d 21 f(X j ) + d 12 f(X j )d 22 f(X j )) + 2u+d 12 f(X j ) (26) 
From (21), (23), (24) and (26) we obtain the distributional limit of the regularized metric in the form 

ds 2 E -dudV + (1 + ^d n f(X j )u+) 2 dX 2 + (1 + ^d 22 f{X^)u+fdY 2 + 

+ ^d 12 f(X j )Af(X j )u 2 + dXdY + 2u+d 12 f(X j )dXdY + (27) 
+ \{d l2 f{Xi)) 2 u\(dX 2 + dY 2 ). 
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Let us sum up the main result of this section in the following diagram 

reg. 



ds 2 = f{x l )5{u)du 2 - dudv + J2(dx*) 2 



?? 



(27) = lim ds 2 



D'-limit 



£^0 



i\2 



ds 2 = f(x l )p £ (u)du 2 - dudv + J2( dxt ) 



ds 2 = -dudV + (2 £(±^4) - 

i=l 

2 

djV e )dudX3 + J2(9 1 x i £ dX^) 2 



The left hand arrow in the above diagram represents the pullback of the distributional metric (5) under the discon- 
tinuous coordinate change (see (11)) 



x l = X l 



i 



(28) 



v = V + f(X>)0(u) + -d t f(Xi)d\f(Xi)u + . 



Note that (28) is the generalization of the Penrose junction conditions (4) to the case of an arbitrary wave profile 
/ (cf. (16) in [13]). Although undefined within Schwartz distribution theory, this transformation can be interpreted 
consistently as the composition of a rcgularization procedure, a smooth transformation of the regularized metric and 
a distributional limit. It is precisely in this sense that the pullback is given by the right hand side of (27) or for short 



ds 2 = -dudV + {5ij + \u + dijf) (5{ + \u+d{f) dX l dX k 



(29) 



For the special case of an impulsive plane wave f(x, y) — x 2 — y 2 this exactly reduces to (3) and in the general vacuum 
case it is equivalent to (2) in [11]. Also, we note that (29) appeared first as (17) in [13]. The regularization procedure 
employed above justifies a posteriori the derivation of (29) by using formal multiplication rules for distributions (cf. 
the remarks following (4)). 



IV. GENERALIZED COORDINATE TRANSFORMATIONS 



In this section we are going to give a precise meaning to the left hand arrow in the above diagram in terms of 
Colombcau's generalized functions. In order to facilitate the following argument we are going to work in the so called 
'simplified' algebra (cf. [10] for a treatment of impulsive pp-waves in the full Colombeau algebra). We note, however, 
that the considerations to follow retain their validity in the 'full' setting as well. Let O be an open subset of R™. With 
£(Q;R n ) = (C^tyR™))^ 00 ' set 

£ M (tt;M. n ) :={(u £ ) £ e£(n;R n ) : VK CC 0, Vm e N 3p e N : 

sup \\D^u £ (x)\\ ^ 0(e- p ) (e^O)} 
xeK 

J\f(Q; R n ) :={{u £ ) £ ef(0;l") : \/K CC 0,Vm € No e N : 

sup HDMu^aOH =O(e 9 )(e^0)} 

xeK 

where D^ m ' denotes the m-th derivative and || . || is any norm on £((R") m ; M"). Then the simplified Colombeau 
algebra £(ft; R") is defined as the differential factor algebra £ M (0; R n )/Af(il; R n ) (see [18-20]). 

Now let ft, ft' be open subsets of R 4 whose w-projection contain the interval [—1, 0]. It is a straightforward consequence 
of Theorem 1 in [10] and the remarks following (14) that (T £ ) £ and (Tc 1 ) e (if defined on Q') are elements of £m{&\ R 4 ) 
and £m(^';R 4 ), respectively. Let us denote their classes in £(£!;R 4 ) (resp. C/(f2';R 4 )) by T and S. We claim that 
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for a suitable choice of ft, 0', S is the inverse in Q of T, which will justify viewing T as a generalized coordinate 
transformation. 

Although composition of Colombeau generalized functions (through componentwise insertion of representatives) is 
not always well defined, a sufficient condition for this operation to yield a unique Colombeau function was given in 
[20], 7.3.1 and 7.4.1, to wit: Let U E (?(Q';R™), V G G(Q;M. n ), then U o V := [{U e o V e ) e ] is a well defined element of 
Q(tl;M. n ) if there exists a representative (U E ) E of [/ such that 

VK cc n 3 if' cc 0' 3r? > s.t. U £ {K) CK'\fe<rj. (30) 

It follows immediately from the remarks following (14) that S satisfies (30). Concerning T, suppose that (u,v,x,y) 
varies in some K CC ft and set (u,V e ,X £ , Y s ) — T e (u, v, x, y). Then 

u s 

x* = xi(X>,u) =*i + \j J dif(xi(X^r))p e (r)drda. 

— S — € 

Since the left hand side and the argument of dif in this equation are bounded by assumption, X\ remains in some 
compact region L in I 2 . We already know that on any such set L (times some compact u- interval) x\ is uniformly 
bounded. Inserting this into 

u us 

V e =v e (V £ ,Xi,u)- I f{x{(X k £1 s)) Pe { S )d S - J J dif(xi(X*,r))xi(X*,r))p e (r)drds 

— e — e — e 

also establishes the desired boundedness property for V e . We conclude that S is indeed the ^-inverse of T (for suitable 
regions Q, O'). Therefore the left hand arrow in the above diagram, i.e. the (generalized) Penrose junction conditions 
(28) turn out to be the macroscopic aspect (or 'distributional shadow') of the (pullback under the) generalized 
coordinate transformation T. 



V. CONCLUSIONS 



In this work we have studied the discontinuous coordinate transformation relating the distributional with the con- 
tinuous form of the impulsive pp-wave metric. It is indeed possible to interpret this change of coordinates as the 
(distributional) limit of a family of smooth transformations which we obtained by a general regularization proceedurc. 
We emphasize the fact that our results are completely regularization independent even within the maximal class 
of regularizations of the Dirac-5, hence "natural" . Note also that we did not have to impose Einstein's equations 
throughout our analysis, hence our results apply to all pp-waves, i.e. even to the non-vacuum case. 
Physically speaking our regularization approach consists in viewing the impulsive wave as a limiting case of a sandwich 
wave with an arbitrarily regularized wave profile. From this point of view the two forms of the impulsive metric arise 
as the (distributional) limits of the sandwich wave in different coordinate systems. It is a well known fact [22] that 
given a one-parameter family of spacetimes one may obtain different limits by introducing a suitable (parameter 
dependent) change of coordinates. This effect is mainly caused by a different identification of points in the manifold. 
More precisely, the transformation in the limit sends finite points to infinity, hence is unbounded. However, in our 
case the limit of the transformation, although discontinuous, remains bounded and the mentioned ambiguities do not 
arise. This is in fact not surprising since our change of coordinates is adjusted to the geometry of the spacetime which 
(in the sense of geodesic completeness) is nonsingular, allowing only for a finite jump of the geodesies. 
Finally, the family of smoothened transformations constitutes a generalized coordinate transformation in the sense 
of Colombeau, whose macroscopic (i.e. distributional) aspect is precisely given by the (generalized) Penrose junction 
conditions (28). 

These results put the formal calculations (using certain ad-hoc multiplication rules for distributions) on solid ground 
and set up a framework which allows to also mathematically identify the two "physically equivalent" descriptions of 
impulsive waves. 
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